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Abstract 

The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical set- 
ting. The leading corrections to this model for very highly charged but dilute ions in thermal 
equilibrium with a weakly coupled, one-component background plasma are explicitly computed, 
and the subleading corrections shown to be negligibly small. This is done using the effective field 
theory methods advocated by Brown and Yaffe. Thus, corrections to nuclear reaction rates that 
such highly charged ions may undergo can be computed precisely. Moreover, their contribution to 
the equation of state can also be computed with precision. Such analytic results for very strong 
coupling are rarely available, and they can serve as benchmarks for testing computer models in 
this limit. 

PACS numbers: 05.20.-y, ll.lO.Wx, 52.25.-b 
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I. INTRODUCTION AND SUMMARY 



Here we describe a plasma configuration whose exact asymptotic solution can be obtained 
in a strong coupling limit. The solution is given by the ion sphere result presented by Salpeter 
plus a simple smaller correction. This is accornglished by using the effective plasma field 



theory methods advocated by Brown and Yaffe P]. In this field-theory language, the old 
Salpeter result corresponds to the tree approximation and our new correction is the one-loop 
term. In usual perturbative expansions, the tree approximation provides the first, lowest- 
order term for weak coupling. Here, on the contrary, the tree approximation provides the 
leading term for strong coupling, with the corrections of higher order in the inverse coupling. 
This is the only example of which we are aware in which the tree approximation yields the 
strong coupling limit. This strongly coupled system is interesting from a theoretical point 
of view and our results can be used to check numerical methods. 

The plasma consists of very dilute "impurity" ions of very high charge ZpC, Zp ^ 1, in 
thermal equilibrium with a classical, one-component "background" plasma of charge ze and 
number density n, at temperature T = 1/13. The background plasma is neutralized in the 
usual way, and it is dilute. We use rationalized electrostatic units and measure temperature 
in energy units so that the background plasma Debye wave number appears as 

K^=^{zefn. (1.1) 

The internal coupling of the background plasma is described by the dimensionless coupling 
parameter 

The assumed weak coupling of the dilute background plasma is conveyed by 

g<^l. (1.3) 

Although the internal coupling of the background plasma to itself is assumed to be very 
weak and the impurity ions are assumed to be so very dilute that their internal interactions 
are also very small, we shall require that the ionic charge Zp is so great that the coupling 
between the impurity ions and the background plasma is very large. To make this condition 
more precise, we define 

Zp = ^. (1.4) 



which is the magnitude of the impurity charge measured in units of the dilute background 
ionic charge. Then the exphcit condition that we require is that 



Since the hmit that we use may appear to be obscure, we pause to clarify it. Even though 
gZp ^ 1, we assume that g is sufficiently small that g'^Zp <^ 1. We may, for example, take 
(7^0 with g"" Zp = const., and a in the interval 1 < a < 2. Then gZp = const. /(^""^ ^ 1 
while g'^Zp = const, g^^" ^ 1. 

Standard methods express the grand canonical partition function in terms of functional 
integrals. Brown and Yaffe do this, introduce an auxiliary electrostatic potential, and 
integrate out the charged particle degrees of freedom to obtain the effective theory. This 
technique will be described in more detail in Sec. |n] below. The saddle point expansion 
of this form for the grand partition function yields a perturbative expansion, with the tree 
approximation providing the lowest-order term. Here, on the contrary, we express the im- 
purity ion number in terms of an effective field theory realized by a functional integral. The 
saddle point of this form of the functional integral involves a classical field solution driven 
by a strong point charge. 

The result for the impurity ion number reads 



Here Np ~ exp{/5yUp} is the number of impurity ions defined by the chemical potential /ip in 
the absence of the background plasma; keeping this chemical potential fixed, the background 
plasma alters this number to be Np. The added ■ ■ ■ stand for corrections to the analytical 
evaluation of the classical action displayed in the Zp^^ and Zp^^ terms of Eq. p.6|) . The sizes 
of these omitted corrections are compared to the exact numerical evaluation of the action 
in Fig. El below. This figure shows that the relative sizes of these terms are small (<C 1) 
in the limit in which we work {gZ ^ 1). The constant C = 0.8499- ■ ■. The final —gZp/3 
term in the exponent is the relatively small one-loop correction. As shown in detail in the 
discussion leading to Eq. ()3.63p below, the error in the result ()1.6p is of the indicated order 
g"^ Zp = g {gZp) and is thus negligible in the limit g <^ I that concerns us. 



> 1 . 



(1.5) 




(1.6) 
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The number correction p.6|) can be used to construct the grand canonical partition 
function Z for the combined system by integrating the generic relation 

for a = p and using the boundary condition that Ai'p — > as — > — cxd. Since depends 
upon the chemical potential /ip only in the factor ~ exp{/5/ip}, this integration gives 

In Z = A/p + Ar(°) . (1.8) 

Here we have identified the constant of integration, the constant that remains when 
vanishes, to be A^*^°\ the number of background plasma particles in the absence of the 
impurity ions. In our limit in which the background plasma is very weakly coupled, N'^^^ ~ 
exp{/3//} is just the number of non-interacting particles of chemical potential /i. 

The equation of state can be found from the well-known relation for a grand canonical 
ensemble with partition function 

(3pV = \nZ. (1.9) 

However, the grand canonical partition function Z is a function of the temperature and 
chemical potentials and, to obtain the equation of state, we must re-express it in terms of 
the observed, physical particle numbers rather than their chemical potentials. 

To do this, we need to express A^^^-* ~ exp{/5/i} in terms of the true number of background 
particles A^, a number that differs from N^^^ because of the presence of the impurity ions. 
There is a significant difference because, although the impurity ions are few in number, they 
are assumed to be extremely highly charged. We again use the general formula ()1.7p . but 
this time to compute N using the solution ()1.8|1 : 

Ar=|^ + Ar(o). (1.10) 
a pfi 

The measured impurity number A^p does depends upon Pfi because it entails the dimen- 
sionless coupling parameter g defined in Eq. p.2|l . For simplicity of exposition, in that 
definition we used a Debye wave number k that was defined in terms of the true background 
density n. Although the distinction is not important for the leading terms that concern us, 
we nevertheless note that the correct wave number that appears in our functional integral 
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formalism involves the 'bare' number density n^'^'' = N^'^^/V, with g ~ ^/nM ~ exp{/3/i/2}, 
and so 

dg 1 



Hence, 



9- (1-11) 



N = N^^) + ^g^^N,. (1.12) 

Using this relation to determine A^^^-* in terms of the physical quantities and Np places 
relationship (jl.9j) of the pressure to the partition function (jl.8p in the proper form of an 
equation of state. To simply bring out the main point, we include here only the leading 
terms, to obtain 

~ (iV - z /^^ff'^' nA T . (1.13) 



10 

Although the fraction of impurity ions in the plasma Np/N may be quite small, there may be 
a significant pressure modification if Zp is very large. Note that the free particle contribution, 
an additional term of A'p, is omitted here since it is not multiplied by the large factor in the 
term that we have retained. 

The number result fll.6p also directly yields the plasma correction to a nuclear fusion 
rate, since 

where Tc is the nuclear reaction rate for a thermal, Maxwell-Boltzmann distribution of the 
initial (1,2) particles in the absence of the background plasma. We use the notation 1 + 2 to 
denote an effective particle that carries the charge {Z_\+Z2)e. This formula was obtained in 
a different guise by DeWitt, Graboske, and Cooper j3[. The relation of the form (jl.l4p that 
we use to previous results is discussed in detail in the Appendix. The formula holds when 
the Coulomb barrier classical turning point of the nuclear reaction is small in comparison 
with the plasma Debye length. This is spelled out in detail in a recent work by Brown, 
Dooling, and Preston ^ who also show that the result ()1.14j) is valid even if the background 
plasma involves quantum corrections. The conditions needed for the formula (jl.l4p to hold 
are also discussed in the work of Brown and Sawyer 0], although sometimes in a rather 
implicit fashion. This work does show, however, that the result ()1.14|1 is valid if KTmax ^ 1, 
where k = (5e^n is the Debye wave number and r^ax is the turning point radius defined 
by Tmax = 2(e^/47rma;^)^/^ where uj = 2nT/his the imaginary time frequency associated 



with the temperature T. It should be remarked that DeWitt, Graboske, and Cooper 
assumed that the nuclear reaction rate formula ()1.14|) held only if the background plasma 
had a classical character, but that the work of Brown, Dooling, and Preston j^] shows that 
it is valid even if the plasma involves quantum effects. 

Our result for the number corrections presents the plasma correction to the fusion 

rate for our special 

3 
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The first line agrees with Salpeter's calculation [1]; the second is new. Again the correction 
can be large. 

We turn now to describe the basis for these results in detail. 
II. REMEMBRANCE OF THINGS PAST 



To begin, we need to review a simple case of the general plasma effective field theory 
formulation presented by Brown and Yaffe P]- First we note that the grand canonical 
partition function for a one-component classical plasma may be expressed as the functional 
integral (which are discussed in detail, for example, in the first chapter of the book by Brown 

), 

/3 
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Here 



A- 



i^exp|-/3i^ 



(2.1) 



(2.2) 



defines the thermal wave length A of the plasma particles of mass m. These particles have 
a chemical potential /i and spin weight (75 so that their density in the free-particle limit is 
given by 



n 



(0) 



^5 



A-3 ef^^ 



(2.3) 



We use rationalized Gaussian units so that, for example, the Coulomb potential appears as 
(j) = e/Airr. We shall be a little cavalier about the uniform, rigid neutralizing background 
that we tacitly assume to be present. We shall explicitly include its effects when needed. 

The validity of the functional integral representation (jzj.lj) is easy to establish. The second 
part in the exponential is written out in a series so as to produce the fugacity expansion 



oo 
n=0 ■ 

[dx] exp 



(2.4) 



a=l 



This Gaussian functional integral can be performed by the functional integration field vari- 
able translation 

n 



a=l 



An r 



Since 



1 



6{r 



,1 I -a), (2.6) 

and the Laplacian can be freely integrated by parts in the quadratic form x(~V^)x, after 
the translation a Gaussian functional integration appears with quadratic form x'(— V^)x' 
with no coupling linear in x'- The original measure [dx] = [dx'] is taken to include factors 
such that this remaining purely Gaussian function integral is simply unity. For pedagogical 
clarity, we make use of the definition (j2.2j) of the thermal wavelength to write the result of 
these manipulations as 



z 



y- 

ra=0 



exp < — /5 



E 



(27r;i)3 

1 " 
5E 



p1 ^ 1 ^ {ze? 
2m 2 ^-^ 47r|ra — r;,| 

,a=l a,b=l 



(2.7) 



This is precisely the familiar fugacity expansion of the classical grand canonical partition 
function. The diagonal sum where a = 6 in the Coulomb potential must be deleted. This 
omission of the infinite self-energy terms is automatic if the dimensional regularization 
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scheme is employed as advocated by Brown and Yaffe Here we shall instead regu- 
late the theory by (at first implicitly) replacing the point source 6{r — r^) with a source 
Sr^t — To) that has a small extent about and (at first implicitly) removing the self energy 
terms, with the limit ^ 5 finally taken in the subtracted theory. 

The derivative of the logarithm of a grand canonical partition function with respect to a 
chemical potential (times (3) gives the particle number conjugate to that chemical potential. 
Thus, if we temporarily add another particle species p of charge Cp = ZpC to the previous 
functional integral, take the described derivative, and then take the limit in which this new 
species is very dilute, we get the desired functional integral representation for the background 
plasma correction to the new species free particle number relation in the presence of plasma 
interactions. 



To express this more precisely, in Eq. (j2.8p Np = gs^ A^^V exp{/3/ip}, where the subscript 
p is used to indicate that these are the properties of the sparsely populated 'impurity' 
ions of charge Cp = ZpC, with V denoting the system volume. So Eq. ()2.8p describes the 
background plasma correction to the free-particle chemical potential - number relationship 
for these p ions immersed in the weakly-coupled, one-component plasma. The original 
chemical potential derivative that leads to this result entailed a volume integral. In virtue 
of the translational invariance of the background plasma, the result is independent of the 
particular value of the spatial coordinate in the electric potential x(r) in the initial factor, 
and this coordinate may be placed at the origin (as we have done), giving the factor e*'^''^^^"^ 
shown. The volume integral then combines to form the total free-particle number A'^"^ that 
appears as a prefactor. We have now subtracted terms from the second exponential, the 
exponential of the action functional of the background plasma, to remove an overall number 
contribution and to include the effect of the rigid neutralizing background. These same 
subtractions must now be made in the normalizing partition function Z that appears in 
the denominator of Eq. ()2.8p . Thus Z is defined by the functional integral of the second 
exponential that appears in Eq. ()2.8p . The effect of the uniform neutralizing rigid background 




(2.8) 
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charge is contained in the term izePx that is subtracted from the exponential exp{izePx}- 
The additional 1 is subtracted from this exponential for convenience. 
To simplify the notation, we write Eq. (|2.8p as simply 

Np = \ J [dx] e-^W , (2.9) 

where the effective action S[x] contains all the terms in both exponents in Eq. ()2.8p . The 
loop expansion is an expansion about the saddle point of the functional integral. At this 
point, the action S[x] is stationary, and thus the field x at this point obeys the classical 
field equation implied by the stationarity of the action. 

The tree approximation is given by the evaluation of S[x] at the classical solution 

X(r) ^ ^0ei(r) , (2.10) 

namely 




+n [e-^-^ciW _ 1 + f3ze<PM] ' (3ZpeS{r) 0ei(r) 

(2.11) 

whose stationary point defines the classical field equation 

-VVci(r) = zen [e-^""*^'^'") - l] + ZpC 5{v) . (2.12) 

This equation defining the classical potential 0ci(r) is of the familiar Debye-Hiickel form, and 
it could have been written down using simple physical reasoning. However, we have placed it 
in the context of a systematic perturbative expansion in which the error of omitted terms can 
be ascertained. In particular, we shall describe the one-loop correction that is automatically 
produced by our formalism. Moreover, we shall prove that higher-order corrections may be 
neglected. Our approach using controlled approximations in which the error is assessed, and 
making precise evaluations of a well defined perturbative expansions in terms of correctly 
identified coupling parameters, differs in spirit from much of the traditional work in plasma 
physics. For example, although previous work has been done by Vieillefosse 12] on the 
solution of the non-linear Debye-Hiickel equation, this work was not done in the context of 
a systematic, controlled approximation. 

9 



The one-loop correction to this first tree approximation is obtained by writing the func- 
tional integration variable as 

X(r) = 20ci(r) + x'(r), (2.13) 

and expanding the total action in Eq. (j2.9p to quadratic order in the fluctuating field x'- 
Since i0ci obeys the classical field equation, there are no linear terms in x' and we have, to 
quadratic order, 

S[x] = 5'[i0ci] 

+ ^ y"(rfr)x'(r) [-V2 + «:2e-'^-^-W] x'(r) , 

(2.14) 



where 



K^ = P{zeYn (2.15) 



is the squared Debye wave number of the mobile ions. The resulting Gaussian functional 
integral produces an infinite dimensional, Fredholm determinant. In this same one- loop 
order, the normalizing partition function Z is given by the same determinant except that it 
is evaluated at 0ci = 0. Hence, to tree plus one-loop order, 

N, = iVf DetV^ [-V^ + k'] {-S[zM} ■ (2-16) 

III. COMPUTATION 
A. Tree 

To solve the classical field equation (j2.12j) in the large Zp limit, we first note that the 
classical potential must vanish asymptotically so as to ensure that the resulting total charge 
density vanishes at large distances form the 'external' point charge Cp = ZpC, 

|r| ^ oo : en [l - e'^^'^'^-'W] ^ . (3.1) 

Since (pd vanishes asymptotically, its defining differential equation ()2.12|1 reduces at large 
distances to 

-\/^^^,{r)c^-K^(j)4v), (3.2) 
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e-K|r| 



and thus, for |r| large, 

0ci(r) ~ (const) ^— — . (3.3) 
|r| 

Since this is exponentially damped, the coordinate integral of the left-hand side of Eq. ()2.12p 
vanishes by Gauss' theorem, and we obtain the integral constraint 



zn 



(dr) [1 _e-/^^^'^^iW] = Zp. (3.4) 



For small r = |r|, the point source driving term in the classical field equation dominates, 
giving the Coulomb potential solution 

0ci(r) ^ ^ . (3.5) 

Thus we write 

<l>M = ^u{0, (3.6) 

where 

e = «:r, (3.7) 
and the point driving charge ZpC is now conveyed in the boundary condition 

m(0) = 1. (3.8) 

The other boundary condition is the previously noted large r limit (|3.3|) which now appears 

as 

^ ^ oo : m(0 ~ e-«. (3.9) 

The action (j2.1H) corresponding to the classical solution is divergent since it includes 
the infinite self-energy of the point charge Cp = ZpC impurity. This self-energy must be 
subtracted to yield the finite, physical action. Following standard practice in quantum 
field theory, the divergent classical action 1)2.111) and the self-energy are first regularized 
— rendered finite — by replacing the point charge with a finite source. The self-energy 
is then subtracted, and finally the point source limit is taken. Regularization is achieved 
by the replacement 6{r) ^^(r), where 5/j(r) is a smooth function of compact support. 
The regularized action obtained by making this substitution in the action ^[i^ci] defined by 
Eq. ()2.11|) will be denoted as S^eg- The regularized self field 0seif(i") is the solution of 

-V2 0seif(r) = V<^i?(r), (3.10) 
11 



and it defines the self-action 

'S'self = 



-/? j (dh) |i(v0scif(r))'- V<^«(r)</>scif(r) 



(3.11) 



Tlie identity 



P j {£r) |(v0scif(r))'-Zpe5H(r)0scif(r)| =0, (3.12) 

whicli is easily verified through partial integration and use of the field equation obeyed by 
0seif can be used to write the self-energy action 1)3.111) as 

Ssei{ = P J (dh) i(^V0,eif(r))' 

= (3 1 {dh) ^ElM, (3.13) 

which is just the impurity's field energy divided by the temperature. It is convenient to use 
this form ()3.13|) in subtracting off the self-energy from S^eg and to also subtract the identity 



(3 J (dh) {V0scif(r) • V0ei(r) - V^«(r) M^)} = 
proved in the same manner as Eq. ()3.13|) . The point source limit S^{r) 
taken to secure the well-defined result 



(3.14) 
6{r) can now be 



S[i4>c\] 
-(3 [ {d'r) 



Ssnh [i(pcl 
1 

2 



VUci(r) 



-n 



where 



^selfV 



Zpe 



(3.15) 
(3.16) 



is the point-source limit of the self-field. 

Using the form ()3.6p for the classical solution we have, remembering that m(0) 

"2 [V(0ei(r)-Cf(r) 
[Zpcf \du{r) 



1 2 



An 



dr 



- - {u{r-) - n(0)) 



du{r-) 
dr 



d 
dr 



1 / \2 

- iu(r) -u(0) 
r V 



(3.17) 
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The final total derivative that appears here gives a null result since the end-point contribu- 
tions vanish. Hence the subtracted action ()3.15|) now appears as 

2 



5'sub[^0cl] 



(3 fdu\ ^ 2 
2 Att \dr J 



exp 



(3ZpZe^ \ ^ (3ZpZe^ 



u > - 



u 



(3.18) 



Changing variables to ^ = nr and using the previously defined plasma coupling constant 

g = /3(ze)^/€/(47r) gives 

72^ / 2 



5'sub[^0cl] 
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exp 



Zp9 



(3.19) 



Requiring that this new form of the action be stationary produces the classical field equation 



exp 



Zpg 



«(0 



Note that the integral constraint ()3.4|) now reads 



oo f 2 

dC- 
9 



1 — exp 



Zp9 



^(0 



(3.20) 



(3.21) 



B. Ion Sphere Model 



In the large Zp limit which concerns us, the short distance form ()3.5|) (multiplied by (3ze) is 
large (compared to one) over a wide range of |r|, and the Boltzmann factor exp{— /?2;e0ci(r)} 
is quite small in this range. We are thus led to the "ion sphere model" brought forth some 
time ago by Salpeter This model makes the step-function approximation 

Zp9 



1 — exp 



(3.22) 



Placing this in the integral constraint 1)3.211) determines the ion sphere radius = /^"^o to 
be given by 

eo = SgZp. (3.23) 
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gZ = 100 




10 



FIG. 1: Numerical solution for u{(^) (solid line), ion sphere model no(C) (short-dashed line), and 
the first correction vi (long-dashed line). For ^ > S^q, uq = ; here = 6.694. 



In the ion sphere model, the classical field equation (j3.2(J|) becomes 



ZpQ 



(3.24) 



and this has the solution, obeying the initial condition Uq{0) = 1, 



MO 



(3.25) 



0, e>eo. 

Here the term linear in ^, a solution of the homogeneous equation, has been determined 
by the continuity at the ion sphere surface, the condition that Mo('Co) = 0. [Without this 
constrain an additional 6{C, — ^o) would appear on the right-hand side of Eq. ()3.24|) .] The 
nature of this "ion-sphere" solution mo(0 together with the exact solution m(^) obtained by 
the numerical integration of Eq. ()3.2()j) . as well as the first correction described below, are 
displayed in Fig. [H 

We have appended the subscript to indicate that this is the solution for the ion sphere 
model. Placing this solution in the new version ()3.19p of the action gives 



-So 



3Z, 
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Zp . 



exp 



1 



(3.26) 

term in the action 



The final —Zp that appears here comes from the 
()3.19|) along with the integral constraint ()3.21|) . This additional —Zp simply adds a constant 
to the chemical potential. Since a constant has no dependence on the thermodynamic 
parameters, this addition has no effect on the equation of state, the internal energy density, 
or any other measurable thermodynamic quantity. Moreover, the contributions of such 
constants clearly cancels in the ratio ()1.14|) that yields the background plasma correction to 
the nuclear reaction rate. 
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C. Ion Sphere Model Corrected 



To find the leading correction to the ion sphere model result, we first cast the exact 
equations in a different form. We start by writing the full solution u{^) as 

u{0 = uo{0 + 4^v{0, (3.27) 

where uo(0 solution ()3.25|) to the ion sphere model ()3.24|) . The exact differential 

equation ()3.20p now reads 



6 







e-^..«o«)/Cexp^-|t;(0[>-^(e-eo) 



(3.28) 



Since Mo(0) = 1 is fixed (reflecting the presence of the large, 'impurity' point charge Zpc), and 
since the solution must vanish at infinity, the proper solution to the non-linear differential 
equation ()3.28j) is defined by the boundary conditions 

i;(0) = 0, e^oo: v{0 ^ ^ ■ (3-29) 

On substituting the decomposition ()3.27p into the action ()3.19|) . the cross term may be 
integrated by parts with no end-point contributions in virtue of the boundary conditions 
(I3.29P on v{^). We take advantage of this to move the derivative of v{^) over to act upon 
uq{^) so that we now have (Puo{^)/d^'^. Using Eq. ()3.24j) for this second derivative and 
identifying the ion sphere part then gives 



^0 f°° 

5'sub[^0ci] = So[i(f)ci] - — d^^v{^) 

Thus far we have made no approximations. To obtain the leading correction to the ion 
sphere result, we note, as we have remarked before, that the factor exp | — ^ "^0(01 is very 
small for ^ < ^o, and so it may be evaluated by expanding mo(0 about ^ = ^o- Using the 
result ()3.25|1 . we find that the leading terms yield 

exp|-^«o(o}^exp|-^ (eo-0'^(eo-o} • (3-31) 
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This approximation is valid for all ^ because when ^ is somewhat smaller than and our 
expansion near the end point breaks down, the argument in the exponent is so large that 
the exponential function essentially vanishes. Indeed, since we consider the limit in which 
.^0 is taken to be very large and the Gaussian contribution is very narrow on the scale set 
by we may approximate 



exp 



Zpg 



uo{i) 



(3.32) 



Here the delta function accounts for the little piece of area that the Gaussian provides near 
the ion sphere radius since 



(3.33) 



With this approximation, an approximation that gives the leading correction for the large 
ZpQ limit in which we work, Eq. ()3.28|) becomes 



-0 (e - eo) f 

^0 



exp 



(3.34) 



It is easy to see that the first correction fi(0 does not alter the integral constraint ()3.21|) . 
Placing the decomposition ()3.27j) in the constraint ()3.2H1 and using the leading-order form 
()3.32|1 together with v{C,) replaced by fi(0 can be used to express the putative change in 
the constraint (j3.2H) in the form 
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exp <J -jvi{^) 



(3.35) 



But Eq. (I3.34|) and partial integration together with the boundary conditions (|3.29|) now 
show that 
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de 



0. 



The — ^o) in Eq- ()3.34|) requires that 



d^ 



d^ 



f=«o-0 



2' 



(3.36) 



(3.37) 
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and 

t^i(eo + 0)-t;i(eo-0)=0. (3.38) 

Since 

and since v{0) =0, we have 

e<eo: v,{0 = c,^, (3.40) 

where ci is a constant that is yet to be determined. For large C,, is small and thus 

obeys the linearized version of Eq. ()3.34j) . 

e»eo: ^%^ = "i(^)' (3.41) 

giving 

e»eo: t7i(0~e-«. (3.42) 

Since this damps rapidly on the scale set by = (3^p5')^^'^ ^ 1, the leading correction f i(0 
that we seek is given by the solution to 

e > eo : = 1 - e--(^) , (3.43) 

which is the previous differential equation ()3.34|) in this region, but with the explicit factors 
of ^/^o and ^q/^ replaced by 1. This new approximate second-order, non-linear differential 
equation is akin to a one- dimensional equation of motion of a particle in a potential with 
C, playing the role of time, and fi(0 playing the role of position. Thus there is an "energy 
constant of the motion". Namely, if we multiply Eq. ()3.43|) by dvi/dC,, we obtain a total 
derivative with respect to ^ whose integral gives 

5(=^)'-M0-e-'« = -l. (3.44) 

where the constant —1 that appears on the right-hand side follows from the limiting form 
as oo. It is easy to show that 

e-'' + v-l>0. (3.45) 
Since asymptotically fi(^) decreases when ^ increases, we must choose the root 
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^2 [e~-^i^) +vi{i)-l\. (3.46) 



The different functional forms for vi{^) in the two regions ^ < and ^ > are joined 
by the continuity constraint ()3.38|) . which we write simply as 

01^0 = ^1(^0), (3.47) 
together with the slope jump ()3.37|) which, using Eq. ()3.46|) . now requires that 

^2 [e-i(5o) + ^;i(eo)-l] = y|e-^^(««) - ^ . (3.48) 

Since we require that ^0 ^ 1, the second term on the right-hand side of this constraint may 
be neglected, which results in a transcendental equation defining ^1(^0), whose solution is 

= 0.6967- ■■ . (3.49) 

We are now in a position to evaluate the leading contribution to the action ()3.3()|1 . Since 
fi(^) damps rapidly on the scale set by ^0, in computing the leading term we can set ^ = ^0 
in the integral that is linear in fi(0- The leading correction is given by 

]-So[iM + Si, (3.50) 

in which 

^2 

5i = -^C, (3.51) 
9 

where 

Here we have omitted the portion 



1 ^?(6 



oj 



. , (3.53) 
^ ^0 

because it is parametrically smaller — it is of relative order l/C,o to the leading terms that 
we retain. We change variables from ^ to vi via 

= dvi , (3.54) 



and use the result ()3.46|1 for the derivative. Hence 



C 



V2 [e-^i - 1] 
1 

2 _ 
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1 rMio) 

+ - / dvi ^2 [e-^'^+Vi - 1] (3.55) 
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1 
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FIG. 2: Ratios of 5[i(/)ci] — Zp for the ion sphere model result (|3.26|) [short-dashed line] and the 
corrected ion sphere model (|3.57() [long-dashed line] to corresponding difference with the action 
(|3.19j) for the exact numerical solution u(^) as functions of gZp. 

is a pure number, 

C = 0.8499---. (3.56) 

In summary, recalling that = (Sfif^p)^^^, we now find that 

(3.57) 

with the leading correction to the ion sphere model exhibited as being of relative order 
1 / {gZp) . Fig. 121 displays the exact numerical evaluation of the action compared with the ion 
sphere approximation [the leading term in Eq. ()3.57|) ] and the corrected ion sphere model 
[the entire Eq. (HTHTj) ]. 




D. One Loop 

The one-loop correction for the background plasma with no "impurity" ions present is 
given by 




(3.58) 



Since we assume that the charge Zp of the "impurity" ions is so large that not only ^ 1, 
but also ZpQ ^ 1 as well, nr^ ^ 1, and the ion sphere radius tq is large in comparison to 
the characteristic distance scale for spatial variation in the background plasma, the Debye 
length . In this case, the term 

exp{-/5ze0(r)} (3.59) 
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in the one-loop determinant that enters into the background plasma correction to the "im- 
purity" number, 

Det-l/2 ^_y2 ^ ^2 g-/3.e</,el] (3 60) 

can be treated as being very slowly varying — essentially a constant — except when it 
appears in a final volume integral. We conclude that in this case of very strong coupling, 

Det^/2 [_V2 + ^2] 




(3.61) 



where in the second equality we have used the ion sphere model that gives the leading term 
for large Z^. 

This result is physically obvious. The impurity ion of very high carves out a hole of 
radius rg in the original, background plasma, a hole that is a vacuum as far as the original 
ions are concerned. The original, background plasma is unchanged outside this hole. This 
ion sphere picture gives the leading terms for very large impurity charge Z^. The corrections 
that smooth out the sharp boundaries in this picture only produce higher-order terms. The 
original, background plasma had a vanishing electrostatic potential everywhere, and the 
potential in the ion sphere picture now vanishes outside the sphere of radius tq. Thus the 
grand potential of the background plasma is now reduced by the amount that was originally 
contained within the sphere of radius ro, and this is exactly what is stated to one-loop order 
in Eq.dSSH). 

This argument carries on to the higher loop terms as well, but we shall now also sketch 
the application of the previous formal manipulations to them as well. 



E. Higher Loops 

As shown in detail in the paper of Brown and Yaffe [^, n-loop terms in the expansion 
of the background plasma partition function with no impurities present involve a factor of 
which combines with other charge and temperature factors to give dimensionless terms 
of the form 

^"-1 I (dh) . (3.62) 
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With the very high Zp impurity ions present, each factor of k is accompanied by 
exp{ — (l/2)/3e 0ci(r)} whose spatial variation can be neglected except in the final, over- 
all volume integral. Thus, in the strong coupling limit of the type that we have set, we have 
the order estimate 

n — loop : 

g^'^K^ j {dh) 1 - exp |-!i±^/5^e0ci(r)| 

g'^-'Khlr^g-Zp. (3.63) 



Again, since we assume that g is sufficiently small so that although gZp ^ 1, g'^ Zp 1, 
all the higher loop terms may be neglected. 

In this discussion, we have glossed over the powers of In g that enter into the higher-order 
terms as well as the quantum corrections that can occur in higher orders. They vanish in 
our strong coupling limit. 
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APPENDIX: RATE RELATED TO PREVIOUS WORK 

n 

We write the result p.l4|) in the form used by Brown, Dooling, and Preston [4] (BDP) 
which is not the notation of DeWitt, Graboske, and Cooper (DGC). In the grand canon- 
ical methods employed by BDP, the temperature and chemical potentials are the basic, fun- 
damental parameters. Thus, in this grand canonical description, the effect of the background 
plasma on nuclear reaction rates appears in terms of number changes with the chemical po- 
tentials held fixed. On the other hand, in the canonical ensemble description employed by 
DGC, the temperature and particle numbers are the basic, fundamental parameters. 

To connect the two approaches, for the relevant case in which "impurity" ions p are 
dilutely mixed in a background plasma, we first note the general structure in the grand 
canonical method. Since the impurities are very dilute, the effect of the background plasma 
on their number is entirely contained in the first term of the fugacity expansion, the linear 
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term in Zp — exp{/9//p}. In the free-particle limit where there is no coupling of the impurities 
to the background plasma, the impurity number density - chemical potential connection 
reads 

<^ = A/ (A.l) 

where Qsj, and Xp are the impurities' spin weight and thermal wavelength, respectively. Thus 
the effect of the background plasma appears as 

Up = e^^ = gs, X;' e^^^ e^- , (A.2) 

where we have chosen to write the plasma correction in terms of an exponential. The only 
feature of the correction Ap that we need note is that it is independent of the impurity 
fugacity Zp since we are working in the — > limit. In summary, the correction in the 
grand canonical description appears as 

=e^^ (A.3) 



(0) 

rVp 



with the total number Np — UpV, where V is the volume of the system. 

The grand canonical partition function Z for the complete system including the various 
impurity ions defines the thermodynamic potential Vt{(5,^) via 

Z = e-^^ , (A.4) 

and the particle number Ng, of species a with chemical potential is given by 

iV. = -|^. (A.5) 

Hence, since generically dN^^'' Idji — PN^^^ , this can be integrated to produce 

where fls is the thermodynamic potential of the background plasma in the absence of the 
extra impurity ions and where, as we have just shown, 

N^'^ = Np exp{-Ap} . (A.7) 

The canonical partition function Zj^ defines the Helmholtz free energy F{(3, N) via 

Zn^c-"^, (A.8) 
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with the connection 

F = n + \2 flaNa. (A.9) 
' 'a 

Since 

/?/i, = ln(nf Aj^?;;) , (A.IO) 
the Helmholtz free energy for a free gas of impurities is thus given by 

iVW) = ivf [in (nW ^;/) - l] . (A.ll) 

The additional ionic impurities change the background plasma free energy from 

Fs = nis + J2^'-Na, (A.12) 

where the sum runs over all the particles in the plasma except for the impurity ions, to 

Using Eq's. (|Xl2)l . (TO)) , and ^K^i produces 

PAFp = Pf^pNp - AtW e^'^ - (3Ff\(3, N,) , (A. 14) 

and, since (ip is fixed in terms of the free gas number densities nf^ = Up exp{— Ap}, we find 
that 

pAFp = -NpAp. (A. 15) 

Thus, in the canonical ensemble approach employed by DGC J^, the previous number ratio 
is expressed in terms of a Helmholtz free energy change, 

""p -exp 1-/3^1 . (A. 16) 



" \ Npj 

These authors sometimes write this in terms of a 'chemical potential'. However, within the 
grand canonical description that we always employ, a chemical potential is an independent 
variable that is not changed as interactions are altered, and so in the context that we use 
this nomenclature is not suitable.. 
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